data {
  int<lower=0> N;
  int<lower=0,upper=1> useforce[N];
  vector[N] igo_pb_use;
  vector[N] salint;
  vector[N] saltan;
  vector[N] terriss;
  vector[N] jointpol7;
  vector[N] rivalry_th;
  vector[N] absidealdiff;
  vector[N] atopally;
}
parameters {
  vector[9] beta;
}
model {

  for (k in 1:9)
        beta[k] ~ normal(0, 10);

  useforce ~ bernoulli_logit(beta[1] + beta[2] * igo_pb_use 
                             + beta[3] * salint
                             + beta[4] * saltan
                             + beta[5] * terriss
                             + beta[6] * jointpol7
                             + beta[7] * rivalry_th
                             + beta[8] * absidealdiff
                             + beta[9] * atopally);
}
generated quantities {
  vector[N] pred;
  for (i in 1:N)
    pred[i] = inv_logit(beta[1] + beta[2] * igo_pb_use[i] 
                             + beta[3] * salint[i]
                             + beta[4] * saltan[i]
                             + beta[5] * terriss[i]
                             + beta[6] * jointpol7[i]
                             + beta[7] * rivalry_th[i]
                             + beta[8] * absidealdiff[i]
                             + beta[9] * atopally[i]);
}